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Abstract 

Onshore  sediment  transport  and  sandbar  migration  are  important  to  the  morphological 
evolution  of  beaches,  but  are  not  understood  well.  Here,  a  new  model  that  accounts  for 
accelerations  of  wave-orbital  velocities  predicts  onshore  sandbar  migration  observed  on 
an  ocean  beach.  In  both  the  observations  and  the  model,  the  location  of  the  maximum 
acceleration-induced  transport  moves  shoreward  with  the  sandbar,  resulting  in  feedback 
between  waves  and  morphology  that  drives  the  bar  shoreward  until  conditions  change. 
A  model  that  combines  the  effects  of  transport  by  waves  and  mean  currents  simulates 
both  onshore  and  offshore  bar  migration  observed  over  a  45-day  period.  A  stochastic 
nonlinear  Boussinesq  model  for  the  evolution  of  waves  in  shallow  water  is  coupled  with  the 
wave- acceleration-driven  sediment  transport  model  to  predict  observed  onshore  sediment 
transport  and  sandbar  migration  given  observations  of  the  offshore  wave  field  and  initial 
bathymetry.  The  Boussinesq-wave  model  has  skill  in  predicting  wave  spectra,  as  well  as 
velocity  and  accceleration  statistics  across  the  surfzone,  but  it  underpredicts  acceleration 
skewness  on  top  of  the  sandbar.  As  a  result,  the  coupled  wave-sediment  transport  model 
underpredicts  sediment  transport,  and  thus  fails  to  move  the  sandbar  onshore.  Although 
the  coupled  wave  and  sediment  model  can  be  tuned  to  yield  skillful  predictions  of  onshore 
sandbar  migration,  in  general,  closer  agreement  between  observed  and  modeled  statistics  of 
the  wave  field  is  essential  for  the  successful  application  of  wave  models  to  predict  sediment 
transport. 

Thesis  Supervisor:  Steve  Elgar 
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Chapter  1 


Introduction 


Sand  bars  are  ubiquitous  morphologic  features  of  sandy  beaches.  They  protect  the  adja¬ 
cent  beach  from  direct  wave  attack  and  are  important  expressions  of  sediment  transport 
patterns  in  the  surfzone,  migrating  back  and  forth  across  the  beach  responding  to  and 
affecting  waves  and  currents,  both  of  which  cause  sediment  transport  and  morphological 
evolution,  including  beach  erosion  and  accretion.  Consequently,  understanding  the  pro¬ 
cesses  that  cause  cross-shore  bar  migration  is  necessary  to  model  and  predict  the  behavior 
of  sandy  beaches  on  time  scales  from  days  to  months. 

During  high  energy  wave  events,  strong  offshore  directed  mean  flows  (undertow)  cause 
rapid  offshore  motion  of  shore-parallel  sand  bars  [Thornton  et  al.,  1996,  Gallagher  et  al., 
1998a],  whereas  during  periods  of  low  energy  waves,  bars  move  slowly  onshore,  often 
becoming  crescentic  [Wright  and  Short,  1984,  Lippmann  and  Holman,  1990].  However, 
the  mechanisms  responsible  for  onshore  sediment  transport  and  bar  migration  are  not 
understood,  and  thus  many  state-of-the-art  models  fail  to  predict  the  observed  beach 
changes.  When  bars  move  onshore  undertow  is  usually  weak,  and  thus  oscillatory  flows 
must  be  the  dominant  hydrodynamic  forcing  [Gallagher  et  al.,  1998a,  Elgar  et  al.,  2001]. 

The  goal  of  this  thesis  is  to  understand  and  model  the  wave-orbital  velocity-driven 
mechanisms  responsible  for  onshore  sediment  transport  and  cross-shore  sandbar  migra¬ 
tion  in  the  surfzone.  The  central  hypothesis  to  be  tested  is  that  fluid  accelerations  asso¬ 
ciated  with  pitched-forward  waves  in  the  surf  zone  result  in  onshore  sediment  transport 
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on  a  natural  beach.  Nearly  continuous  observations  of  waves,  near-bottom  velocity,  and 
bathymetry  collected  along  a  cross-shore  transect  on  a  sandy,  barred  beach  near  Duck, 
NC,  during  the  fall  of  1994  are  used  to  pursue  the  following  specific  objectives: 

•  Implement  and  test  an  acceleration-based  sediment  transport  model  that  explains 
onshore  sediment  transport  and  bar  migration; 

•  Combine  an  existing  energetics-type  cross-shore  sediment  transport  model  [Bagnold, 
1966,  Bowen,  1980,  Bailard,  1981],  known  to  have  skill  in  predicting  offshore  bar  mi¬ 
gration,  with  the  acceleration-based  sediment  transport  model  to  reproduce  observed 
patterns  of  off  and  onshore  sandbar  migration  during  extended  periods  of  time; 

•  Couple  a  stochastic  Boussinesq  wave-shoaling  model  [Herbers  and  Burton,  1997, 
Herbers  et  al.,  2003]  with  the  acceleration-based  sediment  transport  model  to  repro¬ 
duce  observed  onshore  bar  migration  given  offshore  wave  conditions  and  the  initial 
bathymetry. 

1.1  Thesis  Outline 

The  remainder  of  this  chapter  presents  relevant  theoretical  background  for  the  accomplish¬ 
ment  of  the  objectives  outlined  above,  followed  by  a  Data  and  Methods  section.  Results  are 
introduced  in  Chapters  2  and  3  in  the  form  of  two  independent  articles.  Model  calibration 
is  discussed  in  Appendix  A. 

Chapter  2,  “Wave-Acceleration-Induced  Sediment  Transport  and  Sandbar  Migration,” 
an  extension  of  Hoefel  and  Elgar  [2003],  discusses  the  effects  of  wave-acceleration  on 
onshore  sediment  transport  in  the  surfzone,  and  the  implementation  of  an  acceleration- 
based  sediment  transport  model  that  predicts  onshore  sandbar  migration  observed  on  an 
ocean  beach.  The  model  is  forced  with  observed  near-bottom  fluid  accelerations  and  tested 
by  comparison  with  observed  bathymetric  changes.  The  acceleration-based  model  is  then 
combined  with  an  energetics  sediment  transport  formulation  that  reproduces  observed 
sediment  transport  patterns  over  a  45-day  period  during  which  the  sandbar  moved  offshore 
during  storms,  and  onshore  between  storms. 


Chapter  3,  “Coupled  Boussinesq-Wave  and  Sediment  Transport  Model  Applied  to  On¬ 
shore  Sandbar  Migration,”  describes  the  coupling  of  a  stochastic  Boussinesq  wave  model 
[Herbers  and  Burton,  1997,  Herbers  et  al.,  2003]  to  the  acceleration-based  sediment  trans¬ 
port  model  to  predict  observed  onshore  bar  migration  given  observations  of  the  offshore 
wave  field  and  initial  bathymetry. 


1.2  Background 

The  hydrodynamic  processes  that  contribute  to  cross-shore  sediment  transport  in  the  near¬ 
shore  zone  include  wave-orbital  velocities  over  a  range  of  frequencies,  mean  currents,  and 
wave  breaking-induced  turbulence  [Bowen,  1980,  Bailard,  1981,  Roelvink  and  Stive,  1989, 
Wright  et  al.,  1991,  Thornton  et  al.,  1996,  Gallagher  et  al.,  1998a,  Ruessink  et  al.,  1998]. 

During  storms,  breaking-wave-driven  undertow  carries  suspended  sediment  offshore 
[Dally  and  Dean,  1984,  Thornton  et  al.,  1996,  Gallagher  et  al.,  1998a].  Cross-shore  gradi¬ 
ents  in  undertow,  which  is  maximum  just  onshore  of  the  bar  crest  [Thornton  et  al.,  1996, 
Feddersen  et  al.,  1998,  Gallagher  et  al.,  1998a],  result  in  net  transport  of  sediment  from 
onshore  to  offshore  of  the  bar  crest,  causing  the  bar  to  move  offshore.  As  the  bar  moves  off¬ 
shore,  the  locations  of  wave  breaking  (on  the  bar  crest)  and  maximum  undertow  (onshore 
of  the  bar  crest)  also  move  offshore,  resulting  in  feedback  that  promotes  continued  offshore 
bar  migration  as  long  as  there  is  breaking-induced  undertow.  Energetics-type  sediment 
transport  models  [Bowen,  1980,  Bailard,  1981],  originally  developed  to  describe  transport 
under  unidirectional  river  flow  [Bagnold,  1966],  accurately  predict  offshore  migration  of 
sandbars  observed  during  storms  [Stive,  1986,  Russel  et  al.,  1995,  Thornton  et  al.,  1996, 
Gallagher  et  al.,  1998a]. 

For  low  energy  wave  conditions,  when  the  bar  is  observed  to  move  slowly  onshore, 
wave  breaking  on  the  bar  and  corresponding  offshore  mean  currents  are  reduced,  and 
sediment  transport  is  predominantly  caused  by  oscillatory  flows  [Aagaard  et  al.,  1998, 
Gallagher  et  al.,  1998a].  This  transport  is  usually  attributed  to  the  skewed  wave  orbital 
velocities  associated  with  the  sharp  peaks  and  flat  troughs  of  nonlinear  shallow  water  waves 
[Bailard,  1981].  In  the  absence  of  strong  mean  flows,  a  sediment  transport  formulation 
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that  depends  on  an  odd  power  of  wave-orbital  velocity  predicts  more  transport  under 
the  larger  onshore  velocities  associated  with  the  sharply  peaked  wave  crests  than  under 
the  flat,  broad,  wave  troughs,  when  velocities  are  offshore.  However,  velocity-skewness 
based  sediment  transport  models,  such  as  the  energetics  model,  fail  to  predict  onshore  bar 
migration  events  observed  within  the  surfzone  [Roelvink  and  Stive,  1989,  Russel  et  al., 
1995,  Thornton  et  al.,  1996,  Gallagher  et  al.,  1998a],  mostly  owing  to  the  absence  of 
significant  cross-shore  gradients  of  velocity  skewness  near  the  sand  bar. 

As  waves  shoal,  their  shapes  and  orbital  velocities  evolve  from  skewed  profiles  in  inter¬ 
mediate  water  depths,  to  asymmetric,  pitched-forward  shapes  in  the  surfzone  [Elgar  and 
Guza,  1985,  Elgar  et  al.,  1988].  The  sawtooth-like  shape  of  nearly  breaking  and  broken 
waves  is  associated  with  large  fluid  accelerations  and  decelerations  during  the  passage  of 
the  steep  wave  faces,  followed  by  relatively  smaller  decelerations  during  the  passage  of 
the  gently  sloping  rear  of  the  wave,  producing  a  skewed  acceleration  profile.  Field  ob¬ 
servations  [Hanes  and  Huntley,  1986,  Elgar  et  al.,  2001,  Foster  et  al.,  2002],  laboratory 
experiments  [Madsen,  1974,  Cox  et  al.,  1991,  King,  1991,  Sleath,  1999],  and  numerical 
simulations  [Drake  and  Calantoni,  2001,  Hsu  and  Hanes,  in  preparation]  suggest  that  fluid 
accelerations  may  have  a  significant  effect  on  sediment  transport. 

Two-phase  sheet  flow  simulations  [Hsu  and  Hanes,  in  preparation]  corroborate  previ¬ 
ous  field  observations  [Madsen,  1974,  Foster  et  al.,  2002]  that  indicate  severe  bed  failure 
under  the  large  flow  accelerations,  or  horizontal  pressure  gradients,  that  precede  maxi¬ 
mum  onshore  velocities  of  near-broken  waves  in  the  surfzone.  If  sediment  is  mobilized  by 
accelerating  flows,  it  is  expected  that  transport  will  be  larger  when  velocities  are  onshore 
directed  (just  after  strong  accelerations  and  decelerations)  than  during  offshore  directed 
flows,  resulting  in  net  onshore  transport.  Discrete-particle  computer  simulations  of  sheet 
flow  under  unsteady  oscillatory  flows  [Drake  and  Calantoni,  2001]  support  this  hypoth¬ 
esis  and  led  to  the  parameterization  of  acceleration  effects  of  pitched-forward  waves  in 
sediment  transport  as  a  function  of  a  dimensional  form  of  acceleration  skewness  (dis¬ 
cussed  in  Chapter  2).  Effects  of  flow  acceleration  on  bedload  sediment  transport  have  also 
been  accounted  for  in  a  modified,  time-varying  Shields  paramenter  [Nielsen,  2002,  Nielsen 
and  Callaghan,  2003]  used  in  a  Meyer-Peter-type  transport  formulation  to  yield  sediment 


transport  predictions. 

Observations  along  a  cross-shore  transect  of  the  beach  near  Duck,  N.C.,  extending  from 
the  shoreline  to  8-m  water  depth  [Elgar  et  ah,  2001],  show  that  acceleration  skewness  is 
maximum  near  the  bar  crest  for  a  wide  range  of  hydrodynamic  conditions  and  bar  crest 
positions  (Figure  1-1).  The  distribution  of  cross-shore  gradients  of  acceleration  skewness 
is  consistent  with  erosion  offshore  (negative  gradients)  and  accretion  onshore  (positive 
gradients)  of  the  bar  crest  (Figure  1-ld).  Net  onshore  sandbar  migration  would  result 
when  mean  flows  are  weak  ( e.g .  between  days  50  and  60  in  Figure  1-1).  In  contrast,  the 
correlation  between  bar  crest  position  and  location  of  maximum  velocity  skewness  is  low, 
and  cross-shore  gradients  in  velocity  skewness  could  not  account  for  onshore  bar  migration 
(not  shown,  Elgar  et  al.  [2001]). 

Energetics  Sediment  Transport  Model 

Energetics  sediment  transport  models  are  among  the  most  used  in  the  surfzone.  The  orig¬ 
inal  approach  was  developed  by  Bagnold  [1966]  for  unidirectional  flows,  and  later  it  was 
adapted  for  oscillatory  flows  by  Bowen  [1980]  and  Bailard  [1981].  The  model  relates  sedi¬ 
ment  transport  to  fluid  power,  and  it  can  be  written  such  that  it  explicitly  discriminates 
between  mean-flow,  oscillatory-flow,  and  gravity  induced  suspended  load  and  bedload. 
Each  term  is  expressed  as  a  function  of  odd  powers  of  fluid  velocity,  sediment  character¬ 
istics,  and  efficiency  factors  that  represent  the  ratio  of  energy  expended  in  transporting 
sediment  as  bedload  and  suspended  load  to  the  total  rate  of  energy  production  by  the 
stream,  given  by  and  es,  respectively.  The  total  sediment  transport  rate,  Q,  is  given  by 
[Bagnold,  1966] 

Q  —  Qb  +  Qs  —  (^  _  p^g'j  (tan  (f>-  tan/3  +  ( W/u )  —  tan/3)  W  ^  ^ 

where  subscripts  b  and  s  refer  to  bedload  and  suspended  load,  respectively,  ps  and  p  are 
the  sediment  and  fluid  density,  respectively,  <j>  is  the  internal  friction  angle  of  the  sediment, 
/3  is  the  bed  slope,  W  is  the  sediment  fall  velocity,  and  u  is  the  fluid  velocity.  The  last 
term,  w ,  is  the  rate  of  dissipation  of  energy  by  the  fluid  via  bed  friction  and  sediment  flux. 
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Figure  1-1:  Observed  wave  and  near-bottom  cross-shore  velocity  and  acceleration  statistics 
from  Duck,  N.C.  (a)  Significant  wave  height  of  3-hour  long  records  of  sea-surface  eleva¬ 
tion  fluctuations  observed  in  5-m  water  depth  versus  time.  Contours  of  (b)  mean  current 
(negative  values  are  offshore-directed),  (c)  acceleration  skewness,  and  (d)  cross-shore  gra¬ 
dient  of  the  acceleration  skewness  as  a  function  of  cross-shore  location  and  time.  The 
cross-shore  location  of  the  sandbar  crest  is  indicated  by  the  solid  curve  on  each  contour 
plot.  In  panel  d  yellow-red  contours  imply  accretion  (mostly  onshore  of  the  bar  crest), 
and  blue-green  contours  imply  erosion  (offshore  of  the  bar  crest),  consistent  with  onshore 
motion  of  sediment,  especially  when  mean  currents  (panel  b)  are  weak  (e.g.,  between  days 
50  and  60)  (from  Elgar  et  al.  [2001]). 


For  unsteady,  oscillatory  flows  w  is  a  time-varying  quantity  given  by 

w(t)  =  =  P°f\Ht)\3 

where  r  is  the  bed  shear  stress,  Cj  is  the  drag  coefficient,  and  velocity  vector, 

consists  of  u  and  v  components,  parallel  to  the  x  (cross-shore)  and  y  (alongshore)  axes. 
Following  this,  the  time- aver  aged,  cross-shore  volume  sediment  transport  per  unit  width 
per  unit  time  is  given  by  [Bowen,  1980,  Bailard,  1981] 
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Qx  =  Kb{(\u\2u }  +  (|u|2S>  -  ^<|£|3)}  +  +  <M^>  -  ^tanp(\u\5)}  (1.2) 

where  u  and  u  are  the  mean  and  oscillatory  components  of  the  cross-shore  near-bottom 
velocity,  respectively,  tan  f3  is  the  local  cross-shore  beach  slope,  and  angle  brackets  indicate 
time  average.  The  coefficients  K 5  and  Ks  are 


Kb  = 


P  c  tb  . 
( Ps~p)g  7 tan (j)' 


V  -  a  r  _L 
s  {pa-p)g  fW 


Bailard  [1982]  suggested  that  €5  =  0.13  and  es  —  0.01,  similar  to  the  results  of  Bagnold 
[1966]  for  unidirectional  flows.  Church  and  Thornton  [1993],  Thornton  et  al.  [1996],  and 
Gallagher  et  al.  [1998a]  have  used  e 5  =  0.135,  es  =  0.015,  Cf  —  0.003,  and  tan^5>  =  0.63. 

The  cross-shore  transport,  Qx  (eq.  1.2)  is  a  linear  combination  of  bedload  and  sus¬ 
pended  load,  each  driven  by  skewed  oscillatory  velocities,  mean  currents,  and  gravity. 
The  energetics  model  assumes  an  instantaneous  response  of  the  sediment  to  the  flow  and 
yields  a  vertically  integrated  estimate  of  the  transport.  Although  details  of  the  suspended 
sediment  concentration  are  not  considered,  the  model  uses  the  ratio  between  fluid  velocity 
and  sediment  fall  velocity  as  a  surrogate. 


1,3  Data  and  Methods 

Pertinent  information  about  observations  and  data  processing  techniques  are  described 
in  Chapters  2  and  3.  This  section  describes  general  information  or  details  not  discussed 
elsewhere. 

Data  were  collected  as  part  of  the  Duck94  field  experiment  conducted  on  a  micro- 
tidal  («e,  tidal  range  <  1  m)  ocean  beach  located  near  Duck,  on  the  Outer  Banks  of  North 
Carolina,  during  the  fall  of  1994.  The  beach,  roughly  oriented  north-south,  can  present  one 
or  two  sandbars  located  between  30  and  400  m  offshore.  The  mean  foreshore  slope  of  the 
beach  is  approximately  4°  (1:12)  and  the  mean  slope  offshore  of  the  bars  is  approximately 
0.3°  (1:170)  [Lippmann  et  al.,  1993]. 
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Cross-shore  distance  (FRF  coordinates,  m) 

Figure  1-2:  Bathymetric  map  of  the  field  site  located  near  Duck,  N.C.  Filled  black  circles 
represent  colocated  current  meters,  pressure  gages,  and  altimeters,  filled  white  squares  in 
8-m  depth  represent  an  array  of  pressure  gages  to  measure  incident  waves,  and  the  the 
Army  Corps  of  Engineers  Field  Research  Facility  (FRF)  pier  is  to  the  south.  The  map 
is  based  on  26  cross-shore  bathymetric  surveys  conducted  by  the  FRF  on  20  September 
1994.  The  vertical  datum  is  the  National  Geodetic  Vertical  Datum  (NGVD). 

The  data  axe  unique  in  that  they  include  a  wide  range  of  bathymetric  and  hydrody¬ 
namic  conditions  measured  continuously  for  extended  periods  of  time.  Colocated  sonar 
altimeters,  pressure  sensors,  and  bidirectional  electromagnetic  current  meters,  sampled  at 
2  Hz,  were  deployed  on  fixed  frames  spanning  the  inner  and  outer  surfzones,  from  near 
the  shoreline  to  approximately  8-m  water  depth  (Figure  1-2).  All  data  were  stored  in 
3-hour- long  records.  Entire  3-hour  records  for  a  particular  sensor  were  discarded  if  data 
quality  problems  were  detected,  if  the  sensor  was  temporarily  out  of  the  water  owing  to 
tidal  fluctuation,  or  if  the  sensor  was  buried  or  too  close  to  the  seafloor  due  to  bathymetric 
changes.  Current  meters  were  adjusted  vertically  as  the  bathymetry  evolved  to  maintain 
an  approximately  constant  elevation  of  about  50  cm  above  the  seafloor. 

Acceleration  time  series  were  computed  by  differentiating  measured  velocity  time  se¬ 
ries.  When  applicable,  velocity  and  acceleration  time  series  were  Fourier-filtered.  Statis¬ 
tical  moments,  and  other  averaged  quantities,  such  as  sediment  transport  estimates,  were 
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computed  over  3-hour  long  records.  Data  gaps  were  linearly  interpolated  in  time  and 
space. 

Based  on  sediment  transport  estimates  at  each  sensor  location,  cross-shore  gradients 
in  transport  were  computed  for  the  mid-points  between  sensors,  and  then  interpolated 
back  to  match  sensor  locations  following  Gallagher  et  al.  [1998a]. 

The  seafloor  location  under  each  altimeter  (averaged  over  the  approximately  6-cm  di¬ 
ameter  circular  footprint  of  the  altimeter)  was  estimated  for  every  3-hour  period  [Gallagher 
et  al.,  1996].  Spatially  dense  surveys  were  obtained  as  often  as  daily  using  a  GPS  mounted 
on  the  Coastal  Amphibious  Buggy  (CRAB).  These  surveys  were  used  in  conjunction  with 
altimeter  measurements  to  estimate  sandbar  crest  positions. 

Surface  sediment  samples  collected  along  the  surveyed  beach  transect  show  small  tem¬ 
poral  variation  of  cross-shore  distribution  of  grain  sizes  (not  shown)  but  significant  spatial 
variation,  with  poorly  sorted  medium  sand  in  the  swash  and  inner  surf  zones  and  well- 
sorted  fine  sand  seaward  (Figure  l-3a).  Gallagher  et  al.  [1998a]  demonstrated  improved 
energetics  model  performance  by  taking  into  account  cross-shore  grain  size  variations. 
The  same  procedure  is  adopted  here.  Sediment  fall  velocities  were  computed  according  to 
Jimenez  and  Madsen  [2003]  (Figure  l-3b). 

Limited  observations  show  significant  temporal  and  spatial  cross-  and  alongshore  vari¬ 
ations  of  bedform  distribution  during  the  fall  of  1994  and  1997  [Thornton  et  al.,  1998, 
Gallagher  et  al.,  1998b,  2003]  that  include  the  occurrence  of  plane  beds,  wave-orbital 
ripples  (length  scale  0(0.1  m), vertical  scale  0(0.01  m))  ,  and  megaripples  (length  scale 
0(1  m),  vertical  scale  0(0.01  m)).  For  example,  under  small  waves  ( Hs  <  1  m)  and 
weak  currents  (8  Oct  1994)  megaripples  were  observed  in  the  trough  and  seaward  slope 
of  the  bar  crest,  and  small  ripples  were  observed  on  the  bar  crest,  whereas  200  m  to 
the  south,  in  similar  water  depths,  the  seafloor  displayed  a  uniform  distribution  of  small 
ripples,  even  though  the  large  scale  bathymetry  was  alongshore  homogenous  [Thornton 
et  al.,  1998].  Under  moderate  storm  waves  (Hs  ~  2  m)  and  strong  alongshore  currents 
(0.5  to  1  m/s),  Thornton  et  al.  [1998]  observed  wave-orbital  ripples  superimposed  on 
straight-crested  megaripples  in  the  inner  trough  region,  lunate  megaripples  oriented  in 
the  direction  of  the  alongshore  current  on  the  outer  trough,  and  plane  beds  seaward  of 
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Figure  1-3:  (a)  Mean  grain  size  and  (b)  sediment  fall  velocity,  W,  versus  cross-shore 
position  for  samples  collected  in  September  1994. 

the  bar  crest  (no  observations  were  made  on  the  bar  crest).  Following  a  storm  ( Hs  ~  4 
m),  megaripples  were  observed  along  the  axis  of  a  rip  current,  with  plane  beds  farther 
alongshore  [Thornton  et  ah,  1998]. 

Correlations  of  root-mean-square  (RMS)  bed  roughness  with  mobility  number1  and 
Shields  parameter2  estimated  from  measurements  over  a  500  by  700  m  area  during  the 
fall  of  1997  [Gallagher  et  al.,  2003]  suggest  the  transition  from  megaripples  to  plane  beds 
occurs  at  mobility  numbers  ~  150  and  Shields  parameter  ~  0.5.  However,  no  distinction 
of  regimes  can  be  made  below  those  transitions  because  the  correlation  between  RMS 
roughness  and  mobility  number  or  Shields  parameter  is  poor. 


Mobility  number:  ip  =  (^-i)gi>  >  where  U  and  V  are  the  total  (wave  and  current)  instantaneous  cross- 
and  alongshore  velocity  components,  respectively,  s  is  the  specific  gravity,  D  is  the  mean  grain  diameter, 
g  is  the  acceleration  of  gravity,  and  angle  brackets  denote  time  average. 

2Shields  parameter:  0  —  p^s~\)gD  »  where  r  is  the  bed  shear  stress. 
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Chapter  2 

Wave- Acceleration-Induced 
Sediment  Transport  and  Sandbar 
Migration 
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Abstract 

Onshore  sediment  transport  and  sandbar  migration  are  important  to  the  morphological 
evolution  of  beaches,  but  are  not  understood  well.  Here,  a  model  that  accounts  for  fluid 
accelerations  in  waves  predicts  the  onshore  sandbar  migration  observed  on  an  ocean  beach. 
In  both  the  observations  and  the  model,  the  location  of  the  maximum  acceleration-induced 
transport  moves  shoreward  with  the  sandbar,  resulting  in  feedback  between  waves  and 
morphology  that  drives  the  bar  shoreward  until  conditions  change.  A  model  that  combines 
the  effects  of  transport  by  waves  and  mean  currents  simulated  both  onshore  and  offshore 
bar  migration  observed  over  a  45-day  period. 
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2.1  Introduction 


Surfzone  sand  bars  protect  beaches  from  wave  attack,  and  are  a  primary  expression  of 
cross-shore  sediment  transport.  During  storms,  intense  wave  breaking  on  the  bar  crest 
drives  strong  undertow  that  carries  sediment  seaward,  resulting  in  offshore  sandbar  migra¬ 
tion  [Thornton  et  ah,  1996,  Gallagher  et  ah,  1998a]  (Figure  2-la).  If  the  beach  morphology 
is  in  equilibrium,  the  offshore  migration  is  balanced  by  slower  onshore  transport  between 
storms  [Aubrey,  1979,  Wright  and  Short,  1984].  However,  the  causes  of  shoreward  sedi¬ 
ment  transport  and  sandbar  migration  are  not  known,  and  thus  models  for  beach  evolution 
are  not  accurate  [Roelvink  and  Stive,  1989,  Wright  et  al.,  1991,  Thornton  et  al.,  1996,  Gal¬ 
lagher  et  al.,  1998a]. 

As  waves  enter  shallow  water,  their  shapes  evolve  from  sinusoidal  to  peaky,  with  sharp 
wave  crests  separated  by  broad,  flat  wave  troughs.  It  has  been  hypothesized  that  the 
larger  onshore  velocities  under  the  peaked  wave  crests  transport  more  sediment  than  the 
offshore  velocities  under  the  troughs  [Bowen,  1980,  Bailard,  1981].  However,  models  that 
account  for  the  onshore-skewed  velocities  do  not  accurately  predict  onshore  bar  migration 
observed  near  the  shoreline  and  in  the  surfzone  [Roelvink  and  Stive,  1989,  Wright  et  al., 
1991,  Thornton  et  al.,  1996,  Gallagher  et  al.,  1998a],  although  skewed  velocities  may  be 
important  outside  the  surfzone  [Trowbridge  and  Young,  1989].  As  waves  continue  to 
shoal  and  break,  they  evolve  from  profiles  with  sharp  peaks  to  asymmetrical,  pitched- 
forward  shapes  with  steep  front  faces.  Water  rapidly  accelerates  under  the  steep  wave 
front,  producing  high  onshore  velocities,  followed  by  smaller  decelerations  under  the  gently 
sloping  rear  of  the  wave  [Elgar  et  al.,  1988,  2001]  (Figure  2-lb).  Large  accelerations 
generate  strong  horizontal  pressure  gradients  that  act  on  the  sediment  [Madsen,  1974, 
Drake  and  Calantoni,  2001,  Nielsen,  2002].  Although  the  precise  mechanisms  are  not  fully 
understood,  it  has  been  hypothesized  that  if  accelerations  increase  the  amount  of  sediment 
in  motion  [Hallermeier,  1982,  Hanes  and  Huntley,  1986,  Drake  and  Calantoni,  2001,  Elgar 
et  al.,  2001],  there  will  be  more  shoreward  than  seaward  transport  under  pitched-forward 
waves. 
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A  OFFSHORE 

Sandbar  Migration 


Gradients  in  currents  move  the  sandbar  offshore 


Gradients  in  acceleration  skewness  move  the  sandbar  onshore 


Figure  2-1:  Schematic  of  the  feedbacks  that  drive  sandbar  migration,  (a)  Large  waves  in 
storms  break  on  the  sandbar,  driving  a  strong  offshore  directed  current  (undertow)  that  is 
maximum  just  onshore  of  the  bar  crest  [Gallagher  et  al.,  1998a].  The  cross-shore  changes 
(gradients)  in  the  strength  of  the  undertow  result  in  erosion  onshore,  and  deposition 
offshore  of  the  sandbar  crest,  and  thus  offshore  bar  migration.  The  location  of  wave 
breaking  and  the  maximum  of  the  undertow  move  offshore  with  the  sandbar,  resulting  in 
feedback  between  waves,  currents,  and  morphological  change  that  drives  the  bar  offshore 
until  conditions  change,  (b)  Small  waves  do  not  break  on  the  bar,  but  develop  pitched- 
forward  shapes.  Water  rapidly  is  accelerated  toward  the  shore  under  the  steep  front  face  of 
the  waves,  and  decelerates  slowly  under  the  gently  sloping  rear  faces.  Thus,  the  time  series 
of  acceleration  is  skewed,  with  larger  onshore  than  offshore  values  (rectangular  panel).  The 
cross-shore  gradients  in  acceleration  skewness  (maximum  on  the  bar  crest)  result  in  erosion 
offshore,  and  deposition  onshore  of  the  bar  crest,  and  thus  onshore  bar  migration.  The 
location  of  the  peak  in  acceleration  skewness  moves  onshore  with  the  sandbar,  resulting  in 
feedback  between  waves,  currents,  and  morphological  change  that  drives  the  bar  onshore 
until  conditions  change.  23 


2.2  Acceleration-based  Sediment  Transport  Model 


A  surrogate  for  the  effects  of  acceleration  in  pitched-forward  waves  is  a  dimensional  form 
of  acceleration  skewness  [Drake  and  Calantoni,  2001]  ( ie ,  the  difference  in  the  magnitudes 
of  accelerations  under  the  front  and  rear  wave  faces),  a8pike  =  (a3)/{fl2)>  where  a  is 
the  time  series  of  acceleration  and  angle  brackets  denote  averaging.  Discrete-particle 
computer  simulations  of  bedload  transport  driven  by  asymmetrical  waves  characteristic  of 
surfzones  indicate  that  sediment  flux  is  proportional  to  aspike  once  a  threshold  for  sediment 
motion  is  exceeded  [Drake  and  Calantoni,  2001].  Unlike  the  monochromatic  waves  used 
in  the  numerical  simulations,  accelerations  in  random  waves  in  a  natural  surfzone  can 
be  skewed  either  positively  (onshore)  or  negatively  (offshore).  Thus,  the  expression  for 
cross-shore  (x)  acceleration-driven  bedload  sediment  transport  Qa(x)  suggested  by  the 
numerical  simulations  is  extended  to  account  for  random  waves  by  including  a  term  that 
depends  on  the  sign  (ie,  the  direction)  of  aspike >  yielding, 


Qa{%)  ) 


ka  {ci  spike  Sgn  [a,spifce]  CL crit ) 

0 


for  ^  CLcrit 


for  l^spi/cel  ^  CLcrit 


(2.1) 


where  ka  is  a  constant,  sgnfl  is  the  sign  of  the  argument,  and  dcrit  is  a  threshold  that 
must  be  exceeded  for  initiation  of  transport.  By  comparing  model  predictions  with  ob¬ 
servations  (see  Appendix  A  for  discussion  on  model  calibration),  the  optimal  values  of 
ka  =  1.40X1CT4  m  s  and  of  dcrit  —  0.20  m  s'"2  were  determined.  These  parameter  values 
are  within  a  factor  of  5  of  those  suggested  by  the  highly  idealized  discrete  particle  numer¬ 
ical  simulations  [Drake  and  Calantoni,  2001]  ( ka  =  0.26X10  4  ms,  dcrit  =  1-00  ms  ). 
Differences  may  be  owing  to  random  waves,  a  distribution  of  sediment  grain  sizes  and 
shapes,  and  breaking-induced  turbulence  in  the  ocean.  If  it  is  assumed  that  gradients  in 
alongshore  transport  are  negligible,  mass  conservation  in  the  cross-shore  direction  yields 


dh 

dt 


1  dQ( j(x) 


dx 


(2.2) 


where  dh/dt  is  the  change  in  bed  elevation  h  with  time  £,  and  /i  —  0.7  is  a  sediment  packing 
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factor.  Extensions  to  Equation  2.2  to  account  for  alongshore  changes  are  straightforward, 
but  not  necessary  for  the  small  alongshore  gradients  in  transport  inferred  for  the  observa¬ 
tions  discussed  here  [Gallagher  et  ah,  1998a]. 

2,3  Results 

To  test  the  hypothesis  that  the  cross-shore  distribution  of  near-bottom  accelerations  result 
in  overall  onshore  sediment  transport  and  sandbar  migration  when  mean  currents  are  weak, 
morphological  change  predicted  by  the  acceleration-based  model  (Equations  2.1,  2.2)  was 
compared  with  observations  made  along  a  cross-shore  transect  extending  about  400  m 
from  the  shoreline  to  5-m  water  depth  on  the  North  Carolina  coast  [Gallagher  et  al., 
1998a].  The  model  was  initialized  (t  =  0)  with  observed  bathymetry  and  driven  with 
accelerations  observed  with  near-bottom  mounted  current  meters  (Figure  2-2).  During 
a  5-day  period  with  approximately  75-cm-high  waves  and  cross-shore  mean  currents  less 
than  30  cm  s"1,  the  observed  onshore  sandbar  migration  of  about  30  m  was  predicted 
accurately  (Figure  2-2).  A  widely  used  energetics  sediment  transport  model  (Equation 
1.2;  Bagnold  [1966],  Bowen  [1980],  Bailard  [1981],  Thornton  et  ah  [1996],  Gallagher  et  ah 
[1998a])  that  accounts  for  transport  both  by  velocity  skewness  (but  not  acceleration)  and 
by  mean  currents  predicted  no  significant  changes  to  the  cross-shore  depth  profile,  and 
thus  failed  to  predict  the  observed  bar  migration  [Gallagher  et  ah,  1998a]  (not  shown). 
Addition  of  acceleration  effects  (Equation  2.1)  to  the  energetics  sediment  transport  model 
results  in  skilful  predictions  (Figure  2-2). 

During  the  onshore  bar  migration  event,  acceleration  skewness  ( aspike )  increased  from 
small  values  offshore  to  a  maximum  near  the  bar  crest,  and  then  decreased  toward  the 
shoreline,  producing  cross-shore  gradients  in  transport  that  are  consistent  with  erosion 
offshore  and  accretion  onshore  of  the  bar  crest  (Figure  2-3).  The  peak  in  acceleration 
skewness  moved  shoreward  with  the  bar  crest  (Figure  2-3),  resulting  in  feedback  between 
wave  evolution  and  bathymetry  that  promoted  continued  onshore  sediment  transport  and 
bar  movement  until  conditions  changed  (Figure  2-lb).  Feedback  also  occurs  between 
wave-breaking  induced  offshore-directed  mean  currents  (maximum  just  onshore  of  the  bar 
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Figure  2-2:  Observed  and  predicted  cross-shore  bottom  elevation  profiles.  Elevation  of  the 
seafloor  relative  to  mean  sea  level  observed  22  Sep  1994  1900  hrs  EST  (black  solid  curve), 
observed  27  Sep  1994  1900  hrs  (black  dashed  curve),  and  predicted  for  27  Sep  1994  1900  hrs 
by  the  acceleration-based  (red  curve)  and  acceleration  plus  energetics  (using  parameters 
determined  previously  by  Gallagher  et  al.  [1998a],  blue  curve)  sediment  transport  models 
versus  cross-shore  position.  Cross-shore  locations  of  colocated  pressure  sensors,  current 
meters,  and  altimeters  are  indicated  with  triangles,  and  of  colocated  pressure  sensors  and 
current  meters  with  circles.  Observed  near-bottom  velocities  (sampled  at  2  Hz)  were 
low-pass  filtered  (cutoff  frequency  =  0.5  Hz)  and  differentiated  in  time  to  obtain  near¬ 
bottom  acceleration  time  series.  Sediment  transport  fluxes  for  the  model  predictions  were 
computed  from  3-hr  averages  of  observed  near-bottom  velocity  and  acceleration  statistics, 
and  integrated  in  time  with  a  3-hr  time  step  (Equation  2.2)  to  compute  predicted  bottom 
elevation  changes.  Mean  sediment  grain  sizes  ranged  from  0.30  mm  at  the  shoreline  to 
0.15  mm  in  5-m  water  depth  (Figure  1-3). 

crest)  and  morphology  that  results  in  offshore  bar  migration  during  storms  [Thornton 
et  al.,  1996,  Gallagher  et  al.,  1998a]  (Figure  2-la). 

During  an  offshore  sandbar  migration  event  observed  through  a  5-day  storm  (offshore 
wave  heights  >  3  m,  and  undertow  as  strong  as  60  cm  s”1  just  shoreward  of  the  bar  crest), 
the  cross-shore  distribution  of  acceleration  skewness  also  produces  transport  gradients  that 
drive  sediment  from  offshore  to  onshore  of  the  bar  crest  (Figures  1-1  and  2-4).  However, 
during  the  storm  offshore  transport  by  mean  flows  and  skewed  velocities,  predicted  by  the 
energetics  model  (Equation  1.2),  is  larger  than  the  acceleration-induced  transport  (Figure 
2-5),  resulting  in  net  offshore  transport  and  bar  migration  that  is  modeled  accurately  by 
the  energetics  model  (Figure  2-6,  blue  curve).  Addition  of  acceleration- induced  transport 
(Equation  2.1)  to  the  energetics  model  does  not  result  in  significantly  different  predictions 
of  morphological  change  during  the  storm  (Figure  2-6,  red  curve).  However,  during  a 
45-day  observational  period  during  which  the  bar  crest  migrated  offshore  about  130  m 
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Figure  2-3:  Acceleration  skewness  and  bottom  elevation  profiles  during  an  onshore  sandbar 
migration  event,  (a)  Observed  acceleration  skewness  ( aspike )?  (b)  cross-shore  gradient  of 
acceleration  skewness,  and  (c)  seafloor  elevation  relative  to  mean  sea  level  versus  cross¬ 
shore  position.  The  solid  curves  are  observations  from  22  Sep  1900-2200  hrs,  dashed  curves 
are  24  Sep  1300-1600  hrs,  and  dotted  curves  are  27  Sep  1900-2200  hrs. 


during  storms  and  onshore  about  40  m  when  waves  and  mean  flows  were  small  (Figure 
2-7a),  resulting  in  a  net  offshore  migration  of  90  m,  the  combined  energetics  plus  acceler¬ 
ation  model  has  significantly  higher  skill  than  energetics  alone  in  predicting  the  evolution 
of  the  cross-shore  profile.  Although  the  energetics  model  without  acceleration-induced 
transport  predicted  the  offshore  migration,  it  had  limited  skill  predicting  the  total  change 
to  the  beach  over  45  days  because  it  failed  to  predict  onshore  migration  between  storms 
[Gallagher  et  al.,  1998a].  The  energetics  model  extended  to  include  acceleration  effects 
does  better  predicting  the  change  in  the  seafloor  both  onshore  and  offshore  of  the  bar  crest 
(Figure  2-7),  and  the  overall  evolution  of  the  cross-shore  bottom  elevation  profile  (Figure 
2-8). 
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Figure  2-4:  Acceleration  skewness  and  bottom  elevation  profiles  during  an  offshore  sandbar 
migration  event,  (a)  Observed  acceleration  skewness  ( aspike ),  (b)  cross-shore  gradient  of 
acceleration  skewness,  and  (c)  seafloor  elevation  relative  to  mean  sea  level  versus  cross¬ 
shore  position.  The  solid  curves  are  observations  from  10  Oct  2200-0100  hrs,  dashed  curves 
are  14  Oct  0400-0700  hrs,  and  dotted  curves  are  15  Oct  2200-0100  hrs. 

2.4  Discussion 

The  development  of  sediment  transport  models  has  strong  empirical  character  and  re¬ 
lies  on  physical  insights  and  quantitative  data  obtained  in  laboratory  and  field  studies 
[Bowen,  1980,  Bailard,  1981,  Hallermeier,  1982,  Trowbridge  and  Young,  1989,  Dibajnia 
and  Watanabe,  1992,  Thornton  et  al.,  1996,  Gallagher  et  al.,  1998a,  Ribberink,  1998, 
Dohmen- Janssen,  1999,  Malarkey  et  al.,  2003,  Nielsen  and  Callaghan,  2003]. 

Onshore  sediment  transport  under  oscillatory  flows  has  been  attributed  to  skewed 
(sharp-peaks)  near-bottom  wave-orbital  velocities  in  the  shoaling  and  surf  zones  [Bowen, 
1980,  Bailard,  1981,  Dibajnia  and  Watanabe,  1992,  Ribberink,  1998].  Laboratory  mea¬ 
surements  of  transport  of  fine  (0.15  <  D50  <  0.20  mm)  sand  by  second-order  Stokes 
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Figure  2-5:  Total  bottom  elevation  changes  observed  and  predicted  by  the  combined 
energetics  and  acceleration  model  between  10  Oct  1994,  2200  hours  and  15  Oct  1994,  2200 
hours  versus  cross-shore  position.  Observed  bottom  elevation  changes  are  shown  by  the 
dotted  black  curve,  and  changes  predicted  by  individual  sediment  transport  model  terms 
are  shown  by  solid  curves:  mean-flow-driven  changes  (green),  velocity  skewness-driven 
changes  (blue),  acceleration  skewness-driven  changes  (red),  and  total  changes  (black). 

(skewed)  waves  is  consistent  with  velocity  skewness-driven  models,  which  relate  sediment 
transport  to  an  odd  power  of  velocity  [Ribberink,  1998,  Hassan,  2003],  but  transport  of 
fine  (D50  =  0.24  mm)  and  medium  (D50  =  0.44  mm)  sand  under  asymmetrical  (pitched- 
forward)  waves  is  not  [King,  1991,  Ribberink  et  al.,  2000].  For  a  perfectly  asymmetrical 
wave  (ie,  no  skewness),  time  averaged  odd  powers  of  velocity  would  be  zero,  and  thus 
velocity-skewness  based  transport  formulations  would  predict  no  transport.  Sand  trans¬ 
port  under  pitched-forward  waves  is  modeled  more  accurately  by  accounting  for  fluid 
accelerations  [Nielsen,  1992,  2002].  During  the  passage  of  the  steep  front  face  of  a  pitched- 
forward  wave,  abrupt  accelerations  that  precede  maximum  onshore  velocities  produce  thin¬ 
ner  boundary  layers  (ie,  enhanced  bed  shear  stress)  and  greater  pressure  gradients  than 
those  produced  during  the  passage  of  the  gently  sloping  rear  face  of  the  wave.  A  modified 
shear  stress  formulation  that  includes  these  fluid  acceleration  effects  describes  sediment 
transport  in  the  swashzone  [Nielsen,  2002],  where  waves  can  have  strongly  skewed  accelera¬ 
tions  [Raubenheimer,  2002],  and  explains  bedload  transport  rates  of  fine  sand  (D 50  =  0.24 
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Figure  2-6:  Observed  and  predicted  cross-shore  bottom  elevation  profiles.  Seafloor  el¬ 
evation  relative  to  mean  sea  level  observed  10  Oct  1994,  2200  hrs  (solid  black  curve), 
observed  15  Oct  1994,  2200  hrs  (dashed  black),  and  predicted  for  15  Oct  1994,  2200  hrs 
by  the  energetics  (blue)  and  energetics  plus  acceleration  (red)  models  versus  cross-shore 
position. 


mm)  observed  under  pitched- forward  laboratory  waves  [Nielsen  and  Callaghan,  2003]. 

In  the  surfzone,  waves  are  skewed  as  well  as  asymmetric  [Elgar  and  Guza,  1985,  El¬ 
gar,  1987],  including  during  the  onshore  sandbar  migration  observed  between  22  and  27 
September  1994  (Figures  2-2  and  2-9a).  However,  morphological  change  requires  gra¬ 
dients  in  sediment  transport,  and  thus  depends  primarily  on  the  spatial  distribution  of 
hydrodynamic  forcing.  During  the  onshore  sandbar  migration  event,  cross-shore  gradients 
in  normalized  velocity  skewness  were  small  and  uncorrelated  with  sandbar  crest  position 
(Figure  2-9b),  and  therefore  were  unlikely  to  have  caused  net  onshore  sediment  transport 
and  bar  migration.  In  contrast,  gradients  in  normalized  acceleration  skewness  were  larger 
than  those  in  velocity  skewness,  and  imply  onshore  bar  migration  (Figure  2-9b). 

The  precise  mechanisms  of  sediment  transport  that  result  in  the  observed  onshore 
sandbar  migration  are  not  understood.  Observations  by  SCUBA  divers  suggest  there  was 
little  suspended  sediment,  and  none  more  than  a  few  cm  above  the  seafloor,  consistent  with 
measurements  for  a  range  of  conditions  at  this  site  [Conley  and  Beach,  2003].  Immediately 
before  the  sandbar  started  to  migrate  shoreward,  the  seafloor  was  smoothed  by  strong 
currents  (offshore  wave  height  about  2.3  m),  and  remained  relatively  smooth  for  at  least  a 
50-m  wide  cross-shore  section  centered  on  the  bar  crest  during  the  onshore  bar  migration 
(wave  heights  about  0.75  m)  (Figure  2-10).  It  is  hypothesized  that  bedload  dominated 
sediment  transport  when  the  sandbar  migrated  onshore. 
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Figure  2-7:  Observed  wave  height,  cross-shore  sandbar  crest  position,  and  observed  and 
predicted  bottom  elevation  changes  at  four  cross-shore  locations  between  01  Sep  1900  hrs 
and  15  Oct  2200  hrs.  (a)  Significant  wave  height  (4  times  the  standard  deviation  of  3- hr 
long  records  of  sea-surface  elevation  fluctuations  in  the  frequency  bands  between  0.01  and 
0.3  Hz)  observed  in  5-m  water  depth  and  (b)  cross-shore  position  of  the  sandbar  crest 
versus  time.  The  bar  crest  position  was  estimated  from  spatially  dense  surveys  conducted 
with  an  amphibious  vehicle  approximately  bi-weekly,  combined  with  3-hour-long  estimates 
of  seafloor  elevation  from  altimeter  measurements  [Gallagher  et  al.,  1998a]  (Figure  2-2). 
The  shoreline  fluctuated  (owing  to  a  1  m  tide  range)  about  cross-shore  location  £=125 
m.  Observed  (black  circles)  and  predicted  (blue  curve  for  energetics  model,  red  curve  for 
combined  energetics  and  acceleration  model)  cumulative  change  in  seafloor  elevation  at 
cross-shore  locations  (c)  £=161,  (d)  £=220,  (e)  £=265,  and  (f)  £=320  m.  Parameters  in 
the  energetics  model  are  the  same  as  those  in  [Gallagher  et  al.,  1998a]. 
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Figure  2-8:  Observed  and  predicted  cross-shore  bottom  elevation  profiles  spanning  a  45- 
day  period.  Seafloor  elevation  relative  to  mean  sea  level  observed  01  Sep  1900  hrs  (solid 
black  curve),  observed  15  Oct  2200  hrs  (dashed  black),  and  predicted  for  15  Oct  2200  hrs 
by  the  energetics  (blue)  and  energetics  plus  acceleration  (red)  models  versus  cross-shore 
position. 


Although  the  details  of  transport  predicted  by  discrete  particle  simulations  [Drake 
and  Calantoni,  2001]  may  not  be  strictly  applicable  to  fine  grains  transported  in  a  viscous 
regime,  pressure  gradients  produced  by  accelerating  flows  mobilize  and  move  sediment 
regardless  of  grain  size.  Theory  [Madsen,  1974],  numerical  models  [Drake  and  Calantoni, 
2001],  laboratory  measurements  [Sleath,  1999],  and  field  observations  [Foster  et  ah,  2002] 
suggest  that  horizontal  pressure  gradients  can  cause  a  sand  bed  to  become  fluidized  such 
that  resistance  to  stress  is  greatly  reduced,  thus  mobilizing  sediment.  Sediments  in  un¬ 
steady  flows  respond  to  forces  associated  with  flow-induced  drag,  particle  stress  (either 
via  collisions  for  coarse  grains  or  via  vicosity  of  interstitial  fluid  for  concentrated  fine 
grains),  and  horizontal  pressure  gradients  (caused  by  accelerating  flows).  Fine  particles 
respond  rapidly  (order  0.01  s  for  0.2  mm  diameter  sand)  to  both  drag  and  particle  stresses 
in  concentrated  regions.  These  forces  tend  to  balance  each  other,  potentially  allowing 
slower-acting  (order  1  s  for  ocean  waves)  horizontal  pressure  gradient  forces  to  dominate 
[Hsu  and  Hanes,  in  preparation].  In  the  laboratory  [Sleath,  1999]  and  the  field  [Foster 
et  al.,  2002]  blocks  or  plugs  of  fine  sediment  have  been  observed  to  move  under  waves 
with  strong  accelerations.  Under  the  pitched-forward  waves  commonly  observed  in  the 
surfzone,  strong  accelerations  of  near-bottom  fluid  occur  immediately  prior  to  and  during 
onshore-directed  wave  velocities.  Consequently,  mobilized  sediment  will  be  transported 
shoreward,  consistent  with  the  observations. 
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Figure  2-9:  Cross-shore  distribution  of  normalized  fluid  velocity  and  acceleration  skewness 
(calculated  as  the  mean  cube  of  the  demeaned  velocity  and  acceleration  time  series  , 
respectively,  normalized  by  the  variance  [Elgar  and  Guza,  1985]).  (a)  Near-bottom  velocity 
skewness  (dashed  curve)  and  velocity  acceleration  skewness  (solid),  (b)  cross-shore  gradient 
of  velocity  (dashed)  and  acceleration  (solid)  skewness,  and  (c)  elevation  of  the  seafloor 
relative  to  mean  seal  level  versus  cross-shore  position.  Values  are  averages  of  observations 
made  between  22  September  1994  1900  hrs  EST  and  27  September  1994  1900  hrs  ( ie , 
during  the  onshore  sandbar  migration).  Positive  gradients  (panel  b)  imply  accretion  and 
negative  gradients  imply  erosion.  The  vertical  arrows  indicate  the  position  of  the  sandbar 
crest. 
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Figure  2-10:  Roughness  of  the  seafloor  near  the  sandbar.  Root-mean  square  (rms)  seafloor 
roughness  (bedform  heights  are  about  4  times  the  rms  [Gallagher  et  ah,  2003])  versus  time 
before  and  during  the  onshore  bar  migration  (time  =  0  hrs  corresponds  to  22  September 
1994  1900  hrs  EST).  Roughness  was  estimated  as  the  root-mean-square  of  time  series  of 
seafloor  location  determined  by  altimeters  mounted  on  fixed  frames  offshore  (squares), 
on  the  crest  (filled  circles),  onshore  (asterisks),  and  in  the  trough  (open  circles)  of  the 
sandbar  (the  cross-shore  positions  listed  in  the  legend  correspond  to  those  in  Figure  2- 
9).  The  roughness  estimation  assumes  ripples  migrate  beneath  the  altimeters  [Gallagher 
et  ah,  1998b]. 
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Chapter  3 


Coupled  Boussinesq-Wave  and 
Sediment  Transport  Model 
Applied  to  Onshore  Sandbar 
Migration 


Abstract 

A  stochastic  Boussinesq  model  for  the  nonlinear  transformation  of  surface  gravity  waves  is 
coupled  with  a  wave-acceleration-driven  cross-shore  sediment  transport  model  to  predict 
onshore  sediment  transport  and  sandbar  migration  observed  on  a  natural  beach.  The 
coupled  model  is  initialized  with  observed  bathymetry  and  driven  with  wave  frequency- 
directional  spectra  estimated  from  measurements  in  8-m  water  depth.  Similar  to  earlier 
studies,  the  wave  model  has  skill  in  predicting  wave  frequency  spectra,  as  well  as  wave 
velocity  skewness  and  asymmetry  across  the  surfzone.  Dimensional  acceleration  skewness, 
a  third-order  moment  previously  shown  to  explain  onshore  sandbar  migration  in  the  surf 
zone,  also  is  predicted  fairly  well,  although  modeled  values  are  smaller  than  observations 
on  the  crest  of  the  sandbar.  As  a  result,  when  driven  with  modeled  acceleration  skewness, 
the  sediment  transport  model  (calibrated  with  observed  acceleration  moments)  underpre¬ 
dicts  the  rates  and  cross-shore  gradients  of  sediment  transport,  failing  to  move  the  sandbar 
onshore.  Thus,  despite  the  overall  qualitative  agreement  between  observed  and  modeled 
statistics  of  the  wave  field  across  the  surfzone,  closer  agreement  is  essential  for  the  suc¬ 
cessful  application  of  the  Boussinesq-wave  model  to  predict  observed  sediment  transport. 
Improved  skill  of  the  coupled  model  can  be  obtained  by  calibrating  the  sediment  transport 
model  with  modeled  acceleration  moments. 
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3.1  Introduction 


Models  for  the  evolution  of  ocean  surface  gravity  waves  in  shallow  water  are  important 
to  the  prediction  of  nearshore  circulation  and  sediment  transport.  Weakly  nonlinear, 
weakly  dispersive  Boussinesq  wave  transformation  models  accurately  predict  the  observed 
evolution  of  sea-surface  elevation  [Freilich  and  Guza,  1984,  Elgar  and  Guza  1985b,  Elgar 
et  al.,  1990,  Madsen  et  al.  1997,  Norheim  et  al.,  1998,  Herbers  et  al.,  2003,  and  many 
others],  fluid  velocities  [Bosboom  et  al.,  1997],  and  the  frequency-directional  spectrum 
[Norheim  et  al.,  1998,  Herbers  et  al.,  2003]  of  waves  propagating  through  the  shoaling 
and  surf  zones.  In  addition,  Boussinesq  models  have  skill  predicting  third-order  moments 
(skewness  and  asymmetry)  of  velocity  [Elgar  et  al.,  1990,  Bosboom  et  al.,  1997,  Ozanne 
et  al.,  2000,  Herbers  et  al.,  2003]  and  velocity  acceleration  [Elgar  et  al.,  1990]  thought  to 
be  important  to  sediment  transport  [Bowen,  1980,  Bailard,  1981,  Ribberink,  1998,  Drake 
and  Calantoni,  2001,  Elgar  et  al.,  2001,  Hoefel  and  Elgar,  2003]. 

In  shallow  water  (kh  «  1,  where  k  is  the  wavenumber  magnitude  and  h  is  the  water 
depth),  nonzero  third-order  moments  of  wave-orbital  velocity  and  acceleration  arise  from 
near  resonant  triad  nonlinear  interactions  in  which  two  primary  wave  components  with 
frequencies  /i  and  /2  and  wavenumbers  k}  and  fc2  excite  a  secondary  wave  component  with 
the  sum  (h  +/2,  ky+k2  +  6k)  or  difference  (/i  -f2,  kx-k2-  6k)  frequency  and  wavenum¬ 
ber  (where  6k  is  a  slight  mismatch  in  the  wavenumber),  respectively  [Freilich  and  Guza, 
1984],  As  a  result  of  the  energy  transfer  from  the  incident  wave  components  to  higher  and 
lower  frequencies,  the  wave  frequency  spectrum  is  broadened,  and  the  nonlinearly  excited 
secondary  wave  components  are  phase-coupled  to  the  primary  components,  causing  devi¬ 
ations  of  the  wave  field  from  Gaussian  statistics  and  the  steepening  and  pitching-forward 
characteristic  of  near-breaking  waves  [Elgar  and  Guza,  1985a].  Strong  accelerations  un¬ 
der  the  steep  front  of  a  pitched-forward  wave  are  followed  by  smaller  decelerations  under 
the  gently  sloping  rear  of  the  wave,  generating  nonzero  acceleration  skewness  that  has 
been  used  to  explain  wave-driven  onshore  sediment  transport  in  the  surfzone  [Drake  and 
Calantoni,  2001,  Hoefel  and  Elgar,  2003]. 

Here,  a  stochastic  Boussinesq  model  for  the  nonlinear  transformation  of  directionally 
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spread  waves  propagating  over  an  alongshore  uniform  beach  [Herbers  and  Burton,  1997, 
Herbers  et  al.,  2003]  is  coupled  with  a  wave- acceleration-skewness-driven  sediment  trans¬ 
port  model  [Hoefel  and  Elgar,  2003]  to  predict  onshore  sediment  transport  and  sandbar 
migration  observed  on  an  ocean  beach.  When  driven  with  observed  near-bottom  wave- 
orbital  velocities,  the  sediment  transport  model  predicts  the  observed  onshore  sandbar 
migration  [Hoefel  and  Elgar,  2003].  Although  the  Boussinesq  model  qualitatively  predicts 
the  evolution  of  the  wave  field  during  the  5-day  onshore  sandbar  migration  event,  under¬ 
prediction  of  near-bottom  fluid  acceleration  skewness  on  the  sandbar  crest  results  in  poor 
predictions  of  sediment  transport  and  onshore  bar  migration.  The  coupled  model  pre¬ 
dictions  are  improved  by  recalibrating  model  coefficients  to  account  for  the  discrepancies 
between  observed  and  modeled  acceleration  skewness. 

3.2  Boussinesq  Wave  Model 

The  Boussinesq  equations  [Peregrine,  1967]  are  based  on  assumptions  of  weak  nonlinearity 
(a/h  «  1,  where  a  is  the  wave  amplitude)  and  weak  dispersion  (( kh )2  <<  1)  such  that 
the  Ursell  number,  U  —  (a/h)/(k/h)2,  is  approximately  unity.  Stochastic  models  solve 
the  Boussinesq  evolution  equations  for  statistically  averaged  spectral  wave  properties,  and 
thus  are  computationally  efficient.  However,  unlike  deterministic  models  that  solve  the 
approximate  equations  of  motion  without  any  assumptions  about  higher-order  statistics, 
stochastic  models  require  a  statistical  closure  that  may  yield  errors  over  long  propagation 
distances  and  in  regions  of  strong  nonlinearity  [Holloway  and  Hendershott,  1977],  Thus, 
a  relaxation  of  the  bispectrum  to  Gaussian  statistics  over  evolution  distances  comparable 
with  the  surfzone  width  is  used  in  the  closure  to  prevent  spurious  spatial  oscillations  and 
negative  spectral  values  [Orzag,  1970,  Herbers  and  Burton,  1997,  Herbers  et  ah,  2003]. 

Here,  a  two-dimensional  stochastic  Boussinesq  model  for  the  evolution  of  the  frequency- 
directional  spectrum  and  bispectrum  of  surface  gravity  waves  propagating  over  a  gently 
sloping,  alongshore  uniform  beach  [Herbers  and  Burton,  1997]  is  used  to  model  the  evolu¬ 
tion  of  the  wave  field  across  the  shoaling  and  surf  zones.  Dissipation  due  to  wave  breaking 
is  incorporated  with  a  heuristic  frequency-dependent  term  [Kaihatu  and  Kirby,  1995]. 
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The  model,  based  on  a  third-order  closure  (assuming  phase  coupling  between  quar¬ 
tets  of  wave  components  is  weak),  is  cast  in  terms  of  a  coupled  set  of  first-order  evolu¬ 
tion  equations  for  the  frequency-alongshore  wavenumber  (Z)  spectrum  E(fJ)  and  bispec¬ 
trum  J3(/i,Zi,/2,  fe)*  The  two-dimensional  spectrum  !?(/,  Z)  defines  the  energy  density  of 
component  (/,  Z),  and  the  four-dimensional  (complex)  bispectrum  B(f\,  /2,  h)  defines 

the  average  phase  relationship  of  a  triad  consisting  of  components  (/i,/i),  and 

(-/j  -  _  /2).  The  integrals  of  E  and  B  over  all  frequencies  yield,  respectively,  the 

mean  square  and  mean  cube  of  the  sea-surface  elevation.  The  third-order  moments  skew¬ 
ness  and  asymmetry  are  the  integrals  of  the  real  and  imaginary  parts  of  the  bispectrum 
[Elgar  and  Guza,  1985a],  respectively.  Consistent  with  the  depth-averaged  Boussinesq 
equations,  linear  shallow  water  transfer  functions  that  neglect  vertical  variations  in  the 
wave-orbital  velocity  field  were  used  to  convert  spectra  and  bispectra  of  sea-surface  el¬ 
evation  to  spectra  and  bispectra  of  wave-orbital  velocity  and  acceleration,  allowing  the 
third-order  moments  of  velocity  and  acceleration  used  in  the  sediment  transport  model  to 
be  estimated  from  Boussinesq  predictions. 

Discretized  forms  of  the  evolution  equations  are  integrated  from  the  offshore  bound¬ 
ary  across  the  beach  with  a  fixed  step  size  of  2  m.  The  wave  model  is  initialized  with 
bathymetry  observed  at  the  beginning  of  the  onshore  bar  migration  event.  Offshore  bound¬ 
ary  conditions  are  provided  by  wave  frequency-directional  spectra  Eo(f,  l )  estimated  from 
3-hr  long  records  of  near-bottom  pressure  (converted  to  sea-surface  elevation  using  lin¬ 
ear  finite  depth  theory)  measured  with  an  array  of  pressure  gages,  and  the  correspond¬ 
ing  BotfuhJM  approximated  by  second-order  finite-depth  theory  (Hasselmann,  1962, 
Herbers  et  ah,  2003]. 

The  infragravity  frequency  range  /  <  0.05  Hz  was  excluded  from  all  bulk  wave  statistics 
because  strong  shoreline  reflection  of  these  low-frequency  waves  is  not  represented  in  the 
model.  The  maximum  wave  frequency  fmax  was  set  to  0.5  Hz  to  include  higher  frequency 
components  that  may  contribute  to  sediment  transport. 
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3.3  Sediment  Transport  Model 


The  importance  of  wave-acceleration  effects  on  sediment  transport  has  been  suggested 
by  field  observations  [Hanes  and  Huntley,  1986,  Elgar  et  ah,  2001,  Foster  et  al.,  2002], 
laboratory  experiments  [Madsen,  1974,  Cox  et  ah,  1991,  King,  1991,  Sleath,  1999],  and 
numerical  simulations  [Drake  and  Calantoni,  2001,  Hsu  and  Hanes,  in  preparation].  Al¬ 
though  the  precise  mechanisms  involved  are  not  completely  understood,  two-phase  sheet 
flow  simulations  [Hsu  and  Hanes,  in  preparation]  corroborate  previous  field  observations 
[Foster  et  ah,  2002]  that  indicate  severe  bed  failure  under  the  large  flow  accelerations,  or 
horizontal  pressure  gradients,  that  precede  maximum  onshore  velocities  of  pitched- forward 
nearly-broken  or  broken  waves  in  shallow  water.  If  sediment  is  mobilized  by  accelerating 
flows,  it  is  expected  that  transport  will  be  larger  when  velocities  are  onshore  directed  (just 
after  strong  accelerations  and  decelerations)  than  during  offshore  directed  flows,  resulting 
in  net  onshore  transport.  Discrete-particle  computer  simulations  of  sheet  flow  under  un¬ 
steady  oscillatory  flows  [Drake  and  Calantoni,  2001]  support  this  hypothesis  and  led  to  the 
parameterization  of  acceleration  effects  of  pitched-forward  waves  on  sediment  transport 
as  a  function  of  a  dimensional  form  of  acceleration  skewness,  a$Pike  —  (a3) /(a2),  where 
a  is  the  time  series  of  acceleration  and  angle  brackets  denote  averaging.  The  expression 
for  cross-shore  (x)  acceleration-driven  bedload  sediment  transport  Q a  suggested  by  the 
discrete-particle  numerical  simulations  [Drake  and  Calantoni,  2001]  can  be  extended  to 
account  for  random  waves  in  the  surfzone,  yielding  [Hoefel  and  Elgar,  2003] 


Qa{x) 


I^aiOspike  Sgn [dspikel^crit)  for  ^  O crit 

0  for  <  CLcr it 


(3.1) 


where  ka  is  a  constant,  sgn[]  is  the  sign  of  the  argument,  and  a^it  is  a  threshold  that  must 
be  exceeded  for  initiation  of  transport.  By  comparing  model  predictions  with  observations, 
the  optimal  values  of  ka  =  0.014  cm  s  and  of  acru  —  20  cm  s“2  were  determined  [Hoefel  and 
Elgar  2003,  and  also  Appendix  A].  If  it  is  assumed  that  gradients  in  alongshore  transport 
are  negligible  (as  inferred  by  Gallagher  et  al.  [1998]  for  the  observations  discussed  here), 
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mass  conservation  in  the  cross-shore  direction  yields 


dh  _  1  dQa{x) 
dt  [i  dx 


(3.2) 


where  dh/dt  is  the  change  in  bed  elevation  h  with  time  t,  and  fJ.  -  0.7  is  a  sediment 
packing  factor. 


3.4  Field  Observations  and  Data  Analysis 

The  field  observations  discussed  here  were  collected  on  a  barrier  island  beach  in  the  At¬ 
lantic  Ocean,  near  Duck,  N.C.,  during  an  onshore  sandbar  migration  event  observed  be¬ 
tween  22  and  27  September  1994  [Gallagher  et  al.,  1998,  Elgar  et  al.,  2001,  Hoefel  and 
Elgar,  2003].  Offshore  waves  were  approximately  75-cm  high  and  surfzone  cross-shore 
mean  currents  were  less  than  30  cm  s"1.  Twelve  colocated  sonar  altimeters,  pressure 
sensors,  and  bidirectional  electromagnetic  current  meters  (positioned  approximately  50 
cm  above  the  seafloor)  were  mounted  on  fixed  frames  deployed  from  near  the  shoreline 
to  approximately  5-m  water  depth  [Gallagher  et  al.,  1998,  Herbers  et  al.,  2003].  Mean 
sediment  grain  sizes  along  the  transect  ranged  from  0.30  mm  at  the  shoreline  to  0.15  mm 
in  5-m  water  depth  (Figure  1-2,  Gallagher  et  al.  [1998]).  Beach  profiles  used  to  initialize 
the  wave  model  were  obtained  through  linear  interpolation  of  altimeter  measurements. 
Using  high-spatial  resolution  surveys  made  with  an  amphibious  vehicle  does  not  change 
the  results  presented  here. 

All  instruments  were  sampled  at  2  Hz.  Incident  wave  frequency-directional  spectra 
were  estimated  from  3-hr  long  records  from  a  9-element  array  of  bottom-mounted  pres¬ 
sure  gauges  located  in  8-m  water  depth,  about  800  m  from  the  shoreline.  Sea-surface 
elevation  frequency-spectra  across  the  transect  were  estimated  from  3-hr-long  records  of 
near-bottom  pressure  using  linear  finite  depth  theory.  Three-hr  long  records  of  observed 
near-bottom  velocities  were  band-pass  filtered  (with  a  Fourier-type  filter)  between  0.05 
and  0.5  Hz  to  match  wave-model  output,  and  differentiated  in  time  to  obtain  near-bottom 
acceleration  time  series.  Third-order  moments  of  velocity  (skewness  and  asymmetry)  were 
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calculated  as  the  mean  cube  of  the  demeaned  time  series  and  Hilbert-transformed  time 
series,  respectively,  normalized  by  the  variance  [Elgar  and  Guza,  1985a].  Predicted  sed¬ 
iment  transport  fluxes  (Equation  3.1)  at  each  sensor  location  were  computed  from  3-hr 
averages  of  observed  or  modeled  (by  the  Boussinesq  wave  model)  near-bottom  accelera¬ 
tion  statistics,  and  integrated  in  time  with  a  3-hr  time  step  (Equation  3.2)  to  compute 
predicted  bottom  elevation  changes. 

3.5  Coupled  Wave  and  Sediment  Transport  Model  Imple¬ 
mentation 

The  Boussinesq  wave  model  (Section  3.2)  was  run  over  the  initial  observed  bathymetry  for 
the  first  3-hr  boundary  condition  (22  September  1900  hrs).  The  results  were  used  to  update 
the  bathymetry  for  the  next  3-hr  period,  the  same  time  interval  for  which  the  incident 
offshore  wave  spectrum  (estimated  from  observations)  was  updated,  based  on  sediment 
transport  fluxes  (Section  3.3,  Equations  3.1  and  3.2)  computed  from  modeled  aspike  at 
8-m  intervals  along  the  profile  between  the  offshore  pressure-gage  array  (x  =  800  m)  and 
the  sensor  located  closest  to  the  shoreline  (x  =  145  m).  Tests  indicated  no  significant 
difference  between  updating  the  coupled-model  bathymetry  at  spatial  steps  of  8  m  or  at 
the  same  locations  were  sensors  were  deployed. 

The  sequence  of  wave  modeling  and  updating  bathymetry  every  3  hours  was  repeated 
for  the  5  day-long  onshore  migration  event.  The  temporal  evolution  of  the  bathymetry, 
including  the  relatively  slow  onshore  sandbar  migration  is  well  resolved  by  the  3-hr  update 
cycle. 


3.6  Model  Application 

Previous  studies  [Elgar  et.  al,  2001,  Hoefel  and  Elgar,  2003]  suggest  that  onshore  sandbar 
migration  can  result  from  feedback  between  waves  and  morphology.  As  waves  shoal,  veloc¬ 
ity  acceleration  skewness  is  maximum  near  the  crest  of  the  sandbar,  producing  cross-shore 
gradients  in  acceleration-driven  transport  that  result  in  erosion  offshore,  and  accretion  on- 
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shore  of  the  bar  crest,  causing  shoreward  bar  migration.  As  the  sandbar  moves  shoreward, 
so  does  the  location  of  the  maximum  of  acceleration  skewness.  Thus,  accurate  predictions 
of  the  Boussinesq  wave  model  near  the  sandbar  are  of  particular  importance  to  modeling 
morphological  evolution.  Similar  to  previous  studies  [Herbers  and  Burton,  1997,  Herbers 
et  al.,  2003],  the  Boussinesq  wave  model  predicts  sea-surface  elevation  frequency  spectra 
fairly  well  near  the  sandbar  (Figures  3- la  through  3-lc),  except  at  frequencies  above  about 
0.35  Hz.  In  all  3  cases,  the  significant  wave  height  (4  times  the  standard  deviation  of  sea- 
surface  elevation  fluctuations)  was  about  0.9  m  on  the  bar  crest,  which  was  in  about  2-m 
water  depth,  and  thus  there  was  some  wave  breaking.  Deviations  between  modeled  and 
observed  spectral  levels  at  high  frequencies  may  be  owing  to  inaccuracies  in  the  heuristic 
frequency-dependent  dissipation  used  in  the  model,  or  to  neglected  higher-order  nonlinear 
interactions.  The  Boussinesq  model  results  shown  here  (Figures  3-1  through  3-3)  are  from 
the  coupled  wave  and  sediment  transport  model,  and  thus  for  wave  predictions  over  mod¬ 
eled  bathymetry.  However,  discrepancies  between  modeled  and  observed  statistics  of  the 
wave  field  likely  are  not  caused  by  differences  between  modeled  and  observed  bathymetry 
because  Boussinesq  wave  statistics  over  modeled  bathymetry  are  similar  to  Boussinesq 
wave  statistics  over  observed  bathymetry  (not  shown). 

Even  though  the  cross-shore  evolution  of  velocity  skewness  (a  measure  of  the  differ¬ 
ence  in  shape  between  wave  crests  and  troughs,  with  sharp-crested  waves  having  higher 
skewness)  is  reproduced  qualitatively  (Figures  3-ld  through  3-lf)  by  the  Boussinesq  wave 
model,  predicted  values  often  are  higher  than  observed  values,  consistent  with  previous 
findings  ( eg ,  Herbers  et  al.  [2003]).  Velocity  asymmetry  (a  measure  of  the  difference  be¬ 
tween  the  front  and  rear  faces  of  the  waves,  with  pitched-forward  waves  having  higher 
asymmetry)  is  predicted  more  accurately  (Figures  3-lg  through  3-li)  than  velocity  skew¬ 
ness.  Velocity  asymmetry  is  closely  related  to  acceleration  skewness  [Elgar  and  Guza, 
1985a],  but  less  sensitive  to  high  frequency  motions  [Elgar  et  al.,  1990].  Thus,  inaccu¬ 
racies  in  Boussinesq  model  predictions  at  high  frequencies  (eg,  f  >  0.35  Hz  in  Figures 
3-la  through  3-lc)  may  have  a  greater  effect  on  acceleration  skewness  than  on  velocity 
asymmetry  (compare  Figures  3-lg  through  3-li  with  Figures  3-2a  through  3-2c). 

The  cross-shore  distribution  of  dimensional  acceleration  skewness,  aspike >  is  predicted 
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23  Sep 


24  Sep 


27  Sep 


Cross-shore  position  (m)  Cross-shore  position  (m)  Cross-shore  position  (m) 

Figure  3-1:  Observed  (symbols  and  black  curves)  and  modeled  by  the  coupled  Boussinesq 
wave  and  acceleration-induced  sediment  transport  model  (red  curves)  sea-surface  elevation 
spectral  density  near  the  crest  of  the  sandbar  (cross-shore  (x)  positions  are  given  in  each 
panel)  versus  frequency  for  (a)  23  September  1000  hrs,  (b)  24  September  1900  hrs,  and  (c) 
27  September  1000  hrs.  Spectral  estimates  have  approximately  100  degrees  of  freedom  and 
0.0049  Hz  frequency  resolution.  Observed  (black  symbols)  and  modeled  (red  curves)  (d-f) 
near-bottom  velocity  skewness,  (g-i)  near  bottom  velocity  asymmetry,  and  (j-1)  seafloor 
elevation  (relative  to  mean  sea  level)  versus  cross-shore  position  for  the  same  time  periods 
as  the  corresponding  spectra  (a-c). 

fairly  well  by  the  coupled  model  (Figure  3-2a  through  3-2c).  For  example,  on  24  September 
the  agreement  between  observations  and  predictions  is  excellent  (Figure  3-2b).  However, 
often  a spike  is  underpredicted  on  top  of  the  sandbar  crest  (Figures  3-2a  and  3-2c),  some- 
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times  falling  below  the  previously  determined  threshold  for  acceleration- induced  transport, 
cicrit  =  20  cm  s-2  (shaded  area  in  Figures  3-2a  through  3-2c)  [Hoefel  and  Elgar  2003,  and 
also  Appendix  A].  As  a  result,  acceleration-induced  cross-shore  sediment  transport  (Equa¬ 
tion  3.1,  with  ka  —  0.014  cm  s  and  acrit  —  20  cm  s~2)  driven  by  modeled  moments  also  is 
underpredicted  on  the  sandbar  crest  compared  with  transport  driven  by  observed  aspike 
using  the  same  transport  model  coefficients  (Figures  3- 2d  and  3-2f).  Cross-shore  gradi¬ 
ents  of  sediment  transport  estimated  from  modeled  and  observed  aspike  agree  qualitatively, 
and  are  consistent  with  erosion  offshore  (negative  gradients),  and  accretion  onshore  of  the 
sandbar  crest  (positive  gradients,  Figures  3-2g  through  3-2i).  However,  owing  to  the  un¬ 
derprediction  of  aspike  near  the  sandbar  crest,  associated  sediment  transport  gradients 
driven  by  modeled  moments  often  are  smaller  than  those  estimated  from  observations  (eg, 
Figures  3-2g  and  3-2i). 

During  the  5-day  onshore  sandbar  migration  event,  aspike  is  underpredicted  near  the 
sandbar  crest  (Figure  3-3a),  especially  for  the  highest  values  of  aspike ,  which  are  associated 
with  strongly  pitched- forward  waves,  possibly  affected  by  higher-order  nonlinear  interac¬ 
tions  not  included  in  the  Boussinesq  wave  model,  and  by  breaking-induced  dissipation 
that  is  simulated  crudely  in  the  model.  Consequently,  the  coupled  model  often  underpre¬ 
dicts  the  cross-shore  gradients  in  acceleration-induced  sediment  transport  (Figure  3-3b). 
Morphological  evolution  depends  on  gradients  in  sediment  transport  (eg,  Equation  3.2),  so 
that  the  underprediction  of  transport  gradients  by  the  coupled  model  (Figure  3-3b)  results 
in  somewhat  less  erosion  offshore,  and  significantly  less  accretion  onshore  of  the  sandbar 
crest  than  is  observed  and  than  is  predicted  by  gradients  estimated  from  observed  aspike 
(Figure  3-4).  Although  the  coupled  model  predicts  some  of  the  erosion  observed  offshore 
of  the  sandbar  crest  (x  =  265  m,  Figure  3-4b),  it  fails  to  predict  the  accretion  onshore  of 
the  bar  crest  (x  =  220  m,  Figure  3-4a).  As  a  result,  unlike  the  sediment  transport  model 
driven  with  observed  aspikei  the  coupled  model  fails  to  predict  the  onshore  migration  of 
the  sandbar  (Figure  3-5). 

Improved  coupled-model  skill  in  predicting  observed  morphological  change  would  be 
attained  by  improving  the  Boussinesq  model  predictions,  especially  near  the  sandbar  crest, 
perhaps  by  including  higher-order  nonlinearities  or  with  a  better  parameterization  of  wave 
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Figure  3-2:  (a-c)  Observed  (symbols  and  black  curves)  and  modeled  (red  curves)  a  spikes 
associated  (d-f)  acceleration-induced  sediment  transport  Qa  and  (g-i)  cross-shore  gradi¬ 
ents  of  transport  estimated  from  observed  (symbols  and  black  curves)  and  modeled  (red 
curves)  moments,  and  (j-1)  seafloor  elevation  (relative  to  mean  sea  level)  versus  cross¬ 
shore  position.  Left-hand  panels  are  for  23  September  1000  hrs,  middle  panels  are  for  24 
September  1900  hrs,  and  right-hand  panels  are  for  27  September  0100  hrs. 


breaking.  These  improvements  are  beyond  the  scope  of  the  present  study.  Alternatively, 
the  effect  of  underpredicting  aspike  on  the  sandbar  crest  ( eg ,  Figure  3-3a)  can  be  compen¬ 
sated  partially  by  lowering  the  threshold  for  initiation  of  acceleration- induced  transport, 
dcrit  •  Comparison  of  coupled-model  predictions  with  observations  of  morphological  change 
suggest  that  for  the  onshore  migration  event  discussed  here,  the  highest  predictive  skill  is 
attained  for  acrit  =  8  cm  s-2  and  ka  —  0.010  cm  s  (Figure  3-5). 
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Normalized  cross-shore  position  (x/xbar) 


Figure  3-3:  (a)  Dimensional  acceleration  skewness  ( aspike )  and  (b)  cross-shore  gradients 
in  acceleration-driven  sediment  transport  Qa  versus  cross-shore  position  ( x )  normalized 
by  the  position  of  the  sandbar  crest  (xbar)-  Thus,  x/xbar  =  1  is  the  cross-shore  position 
corresponding  to  the  (moving)  bar  crest.  Symbols  are  observed  values,  and  the  red-shaded 
areas  are  the  envelopes  of  values  predicted  by  the  coupled  model  from  41  consecutive  3-hr 
runs  beginning  22  September  1900  hrs  and  ending  27  September  1900  hrs. 
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Figure  3-4:  Observed  (black  circles)  and  predicted  (blue  circles  for  sediment  transport 
model  driven  with  observed  aspike  and  red  circles  for  the  coupled  wave  and  sediment 
transport  model)  cumulative  change  in  seafloor  elevation  versus  time  at  cross-shore  loca¬ 
tions  (a)  x  -220  m,  and  (b)  x  =265  m.  The  dotted  vertical  lines  indicate  the  time  periods 
depicted  in  Figures  3-1  and  3-2:  23  September  1000  hrs,  24  September  1900  hrs,  and  27 
September  0100  hrs. 
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Figure  3-5:  Elevation  of  the  seafloor  (relative  to  mean  sea  level)  versus  cross-shore  position 
observed  on  22  September  1900  hrs  (black  solid  curve),  observed  on  27  September  1900 
hrs  (black  dashed  curve),  predicted  on  27  September  1900  hrs  by  the  acceleration-based 
sediment  transport  model  driven  with  observed  aspike  (blue  solid  curve,  ~  20  cm  s  2, 
ka  -  0.014  cm  s),  predicted  on  27  September  1900  hrs  by  the  coupled  wave  and  sediment 
transport  model  (red  solid  curve,  0,^0,  —  20  cm  s“2,  ka  =  0.014  cm  s),  and  predicted  on 
27  September  1900  hrs  by  the  coupled  wave  and  sediment  transport  model  (red  dashed 
curve,  dcrit  =  8  cm  s"2,  ka  —  0.010  cm  s). 

The  sensitivity  of  the  coupled  wave  and  sediment  transport  model  to  the  choice  of  high 
frequency  cut-off  was  briefly  investigated.  Initially,  the  same  high  frequency  cut-off  (0.5 
Hz)  used  for  the  data-driven  sediment  transport  model  (Chapter  2)  was  chosen  for  the 
Boussinesq  wave  model.  The  performance  of  the  sediment  transport  model  driven  with 
observations  is  insensitive  to  the  choice  of  high  frequency  cut-off  to  the  extent  that  equally 
skillfull  predictions  can  be  obtained  for  different  high  frequency  cut-offs  provided  that 
recalibrated  sediment  transport  model  constants  are  used.  However,  the  underprediction 
of  wave  spectral  densities  by  the  wave  model  at  frequencies  higher  than  about  0.35  Hz 
(Figure  3-la  through  3-lc)  suggests  that  the  wave  model  might  be  sensitive  to  the  choice 
of  high  frequency  cut-off.  Better  agreement  between  observed  and  modeled  moments  is 
obtained  by  lowering  the  high  frequency  cut-off  to  0.4  Hz.  For  this  frequency  cut-off,  both 
the  observation-driven  and  the  coupled  Boussinesq- wave  sediment  transport  models  yield 
reasonable  agreement  with  observed  profile  evolution  using  the  same  sediment  transport 
model  constants  ( acru  =  11  cm  s~2,  ka  =  0.017  cm  s;  Figure  3-6). 
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Figure  3-6:  Elevation  of  the  seafloor  (relative  to  mean  sea  level)  versus  cross-shore  position 
observed  on  22  September  1900  hrs  (black  solid  curve),  observed  on  27  September  1900 
hrs  (black  dashed  curve),  predicted  on  27  September  1900  hrs  by  the  acceleration-based 
sediment  transport  model  driven  with  observed  aspike  (blue  solid  curve),  and  predicted 
on  27  September  1900  hrs  by  the  coupled  wave  and  sediment  transport  model  (red  solid 
curve).  For  both  the  data-driven  and  the  coupled  wave  and  sediment  transport  models 
the  high  frequency  cut-off  was  set  to  0.4  Hz,  and  =  11  cm  s"2  and  ka  =  0.017  cm  s. 

3.7  Summary 

Although  the  cross-shore  evolution  of  waves  shoaling  and  breaking  on  a  barred  beach  is 
modeled  well  by  a  stochastic  Boussinesq  wave  model  initialized  with  the  wave  field  ob¬ 
served  at  the  seaward  edge  of  the  beach  profile  [Herbers  and  Burton,  1997,  Herbers  et  ah, 
2003],  relatively  small  errors  in  predicted  higher-order  moments  of  the  wave  field  have  im¬ 
portant  effects  on  subsequent  sediment  transport  and  morphological  change  modeling.  In 
particular,  underprediction  of  near-bottom  fluid  acceleration  skewness  near  the  crest  of  the 
sandbar  results  in  reduced  skill  in  acceleration-based  sediment  transport  model  predictions 
relative  to  predictions  based  on  observed  acceleration  time  series.  Consequently,  the  cou¬ 
pled  Boussinesq- wave  and  acceleration- induced  sediment  transport  model  introduced  here 
underpredicts  the  onshore  migration  of  a  sandbar  observed  on  an  ocean  beach  for  moder¬ 
ate  wave  conditions  with  weak  mean  cross-shore  currents.  However,  adjustment  of  model 
parameters,  in  particular  lowering  the  threshold  for  initiation  of  acceleration-induced  sed¬ 
iment  transport  relative  to  the  optimal  threshold  determined  from  data-driven  transport 
modeling,  results  in  improved  predictive  skill  of  morphological  change.  Further  research 
is  needed  to  improve  the  predictions  of  the  Boussinesq  wave  model,  especially  in  regions 
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of  strong  nonlinearity  and  wave-breaking  induced  dissipation  (eg,  near  the  sandbar  crest), 
and  to  calibrate  sediment  transport  model  coefficients  with  additional  data  sets.  How¬ 
ever,  the  results  presented  here  suggest  that  coupling  a  stochastic  Boussinesq  wave  model 
with  a  wave-induced  sediment  transport  model  can  result  in  skillful  predictions  of  the  on¬ 
shore  sediment  transport  and  sandbar  migration  observed  on  an  ocean  beach  when  mean 
currents  are  weak. 
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Appendix  A 


Combined  Energetics  and 
Acceleration  Model  Performance 


Model  performance  can  be  quantified  by  means  of  the  normalized  root  mean  square  (RMS) 
prediction  error  Erms  / Arms,  where  Erms  is  the  RMS  error  between  final  predicted  and 
observed  3-hour  averaged  profiles  and  Arms  Is  the  RMS  change  between  the  initial  and 
final  observed  profiles.  Model  skill  is  defined  as  (after  Davis  [1976]  and  Gallagher  et  ah 
[1998]) 


skill  =  1  - 

Arms 


(A.1) 


where, 


Erms  = 


\ 


1  N  2 
iyF  ^  >  \Sn(hpredicted,n  ^ observed, n)] 


N 


71=  1 


(A.2) 


and 


Arms  ~ 


77 


71=1 


(A.3) 


where  N  is  the  number  of  sensors  along  the  cross-shore  transect,  hpre<ncted,n  is  the  predicted 
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final  profile  at  each  location  n,  h^served,n  is  the  observed  final  profile,  and  An  is  the 
measured  change  in  seafloor  elevation  at  each  of  the  N  sensor  locations.  To  account  for 
irregular  sensor  spacing,  estimates  at  each  cross-shore  location  are  weighted  by  Sn  given 
by 


f  ^  (dxn— 1 4~<jxn) 


V  dx 


for  n  =  2,  ...,N  —  1 
for  n=l,N 


(A.4) 


where  dxn  =  xn+i  —  xn  is  the  separation  between  adjacent  sensors  located  at  cross-shore 
positions  xn+ 1  and  xn ,  and  dx  is  the  average  sensor  separation  along  a  transect. 

Model  skill  is  maximum  and  equal  to  1  when  there  is  perfect  agreement  between  final 
observed  and  predicted  profiles.  If  skill  —  0,  prediction  errors  are  as  large  as  observed 
changes,  whereas  if  skill  <  0  prediction  errors  are  larger  than  observed  changes.  If  ob¬ 
served  changes  between  initial  and  final  profiles  are  small  ( ie ,  Arms  is  small),  the  corre¬ 
sponding  model  skill  may  be  low  even  if  errors  (Erms)  in  predictions  are  small  [Gallagher 
et  ah,  1998]. 

The  skill  of  the  combined  acceleration  and  energetics  model  was  evaluated  as  a  function 
of  the  acceleration-based  model  constants  ka  and  dcm  for  three  simulation  periods  for  the 
Duck94  data  set  (Figure  A-l).  These  periods  include  the  onshore  bar  migration  event 
between  22  and  27  Sep  (Figure  2-2),  an  offshore  bar  migration  event  observed  during  a 
storm  between  10  Oct  and  15  Oct  (Figure  2-6),  and  a  45-day  long  simulation  between 
01  Sep  and  15  Oct  that  includes  on  and  offshore  bar  migration  events,  as  well  as  lower 
energy  wave  events  during  which  no  significant  bathymetric  changes  were  observed  (Figure 
2-7).  Calibration  and  model  performance  of  the  energetics  formulation  has  been  evaluated 
previously  by  Thornton  et  al.  [1996],  and  Gallagher  et  al.  [1998]. 

Simulations  for  the  onshore  bar  migration  event  (22  to  27  Sep,  Figure  A-la)  indicate 
little  model  sensitivity  to  ka  (ie,  the  maximum  skill  of  approximately  0.7  is  attained 
for  0.02  <ka  <  0.2  cm  s)  and  large  sensitivity  to  the  acceleration  threshold  for  sediment 
transport,  dcrit  (ie,  maximum  model  skill  is  attained  for  a  narrow  range  of  dent).  A  similar 
sensitivity  to  ka  and  is  observed  for  the  offshore  bar  migration  event  (10  to  15  Oct, 
Figure  A- lb).  Including  the  effects  of  acceleration- induced  transport  (ie,  when  dcrit  <  50 
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(A)  Sep  22  to  Sep  27  (B)  Oct  1 0  to  Oct  1 5  (C)  Sep  01  to  Oct  1 5 


0  20  40  60  80  100  0  20  40  60  80  100  0  20  40  60  80  100 

Acceleration  threshold  (cm/s2)  Acceleration  threshold  (cm/s2)  Acceleration  threshold  (cm/s2) 


Figure  A-l:  Combined  energetics  and  acceleration  sediment  transport  model  skill  as  a 
function  of  ka  and  acrit  for  three  simulations,  (a)  22  Sep  to  27  Sep,  1994  (onshore  bar 
migration  event);  (b)10  Oct  to  15  Oct,  1994  (offshore  bar  migration  event);  and  (c)  Sep 
01  to  Oct  15  (onshore  and  offshore  bar  migration  events).  All  skills  <  0  were  grouped  into 
the  skill  =  0  bin  for  better  visualization. 

cm/s2)  results  in  improved  skill  relative  to  predictions  of  the  energetics  model  alone  (ze, 
when  acrit  >  50  cm/s2).  Simulations  that  include  on  and  offshore  bar  migration  events, 
as  well  as  periods  with  no  significant  bathymetric  changes,  indicate  a  limited  optimum 
parameter  space  centered  around  the  values  of  ka  ~  0.01  cm  s  and  acrit  ~  20  cm/s2  (Figure 
A-lc). 
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